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1.  INTRODUCTION 


This  report  presents  an  analysis  of  selected  long-wavelength  infrared 
(LWIR)  limb  radiance  data  measured  in  the  1983  ELC-1  rocket  experiment. 
From  the  CO2  ^2  band  emission,  information  on  the  mesospheric  temperature 
profile  and  the  radiometer  absolute  calibration  is  extracted.  The  results  are 
combined  with  a  measurement  of  water  vapor  rotational  emission  in  the 
23-29  /xm  region  to  derive  a  water  vapor  profile  for  the  altitude  range  55  km 
-  85  km.  Over  most  of  these  altitudes  the  emission  is  optically  thick.  The 
methods  presented  here  can  also  be  applied  to  data  obtained  in  the 
CIRRIS-IA  experiment  flown  on  the  Space  Shuttle  STS-39  in  April,  1991,(1) 
which  viewed  the  earth  limb  with  near-global  coverage. 

We  begin  with  a  brief  description  of  the  experiment.  In  the  next 

section,  the  temperature  profile  is  discussed.  In  the  subsequent  sections,  we 
discuss  the  data  reduction  and  inversion  methods  used  to  derive  the  water 
vapor  profile.  Finally,  the  ELC-1  water  vapor  profile  is  compared  to 

ground-based  microwave  measurements  and  to  infrared  absorption  measurements 
from  spacecraft. 

The  ELC-1  rocket  experiment  was  launched  on  October  25,  1983  from 

the  White  Sands  Missile  Range,  at  local  time  0426  hours.  The  geomagnetic 

activity  index  Kp  for  that  day  was  1.61  and  the  solar  flux  index  F10.7  was  89. 
The  atmospheric  model  MSIS-90(2)  was  run  with  these  inputs  to  define  the 
profiles  used  in  the  data  analysis. 

The  ELC-1  radiometer  contained  an  array  of  cryogenically-cooled 
detectors  covering  the  spectral  range  4  ^m  -  29  jtim.  The  basic  sensor  was 
previously  flown  in  the  Zodiacal  Infrared  Project  (ZIP).(3)  For  the  ELC-1 
experiment  the  sensor  was  mated  to  a  high-off-axis-rejection  telescope,  and  the 
detectors  were  fitted  with  bandpass  filters  to  yield  the  spectral  responses  given 
in  Table  1. 

The  rocket  reached  an  apogee  of  300  km  at  296  sec  time  after  launch 
(TAL).  Limb-viewing  horizontal,  vertical,  drop  and  roll  scans  were  performed 
during  the  interval  110  -  478  sec  TAL.  The  aspect  solutionC^)  was  obtained 
by  combining  star  camera  and  gyroscope  data  and  iterating  to  obtain 
self-consistent  IR  radiance  profiles.  In  this  report,  data  from  two  sets  of 


scans  are  analyzed:  a  set  of  four  horizontal  scans  taken  during  the  interval 
230  -  256  sec  TAL  and  a  set  of  four  vertical  scans  taken  during  353  - 
478  sec  TAL.  Based  on  the  self-consistency  of  the  radiance  profiles  and  the 
use  of  the  star  data  to  establish  the  absolute  aspect,  we  estimate  the  tangent 
heights  to  be  accurate  to  within  2  km  or  better. 

The  CO2  t'2  band  emissions  were  measured  by  detectors  7,  9  and  31, 
and  the  H2O  pure  rotational  emission  was  measured  by  detectors  27  and  28. 
Other  detectors  provide  useful  data  on  the  O3  1/3  cold  and  hot  bands  in  the 
9-12  nm  region  and  on  the  NO  fundamental  at  5.3  nm.  These  data  will  be 
the  subject  of  a  future  report. 


Table  1.  Bandpasses  for  ELC-1  Detectors 
for  a  Blackbody  at  300K. 


Detector 

Number 

Band  Center 
Wavelength 
(Min) 

Optical 

Passband 

(Mm) 

1,2 

15.137 

14.76-15.52 

3,4 

12.498 

12.02-12.97 

5,6 

9.630 

9.28-9.98 

7,8 

16.435 

15.79-17.08 

9,10 

13.747 

13.06-14.43 

11,12 

9.586 

9.24-9.93 

13,14 

10.583 

10.02-11.15 

15,16 

15.180 

14.82-15.54 

17,18 

14.097 

13.56-14.63 

21,22 

5.571 

5.11-6.03 

23,24 

4.653 

4.34-4.97 

27,28 

26.201 

23.38-29.02 

29,30 

11.499 

11.00-12.00 

31,32 

15.197 

14.81-15.58 

2.  TEMPERATURE  PROFILE 


The  inferred  H2O  profile  is  sensitive  to  the  assumed  temperature 
profile.  Our  approach  for  determining  the  temperature  profile  of  the  ELC-1 
atmosphere  was  a  simple  one.  The  initial  density  and  temperature  profiles 
were  taken  from  the  empirical  model  MSIS-90.(2)  For  the  altitude  range  of 
interest  (up  to  85  km),  the  CO2  mixing  ratio  was  set  to  the  value  at  ground 
level.  The  CO2  spectral  radiance  profile  in  the  u2  region  was  then  calculated 
using  the  Strategic  High  Altitude  Radiance  Code  (SHARC).(5) 

The  assumed  CO2  density  profile  should  be  reasonably  accurate,  and  in 
any  case  the  radiance  is  insensitive  to  it  since  the  CO2  bands  are  optically 
thick.  Therefore,  the  only  significant  sources  of  disagreement  between  the 
SHARC  calculation  and  the  ELC-1  data  would  be  errors  in  the  temperature 
profile  or  in  the  radiometric  calibration.  A  calibration  error  would  cause  the 
radiance  profile  to  be  uniformly  scaled.  Errors  in  the  temperature  profile,  on 
the  other  hand,  could  cause  discrepancies  that  vary  with  tangent  height. 

The  bandpass  radiances  were  calculated  by  SHARC  for  detectors  7,  9, 
and  31.  Detectors  7  and  9  measured  emissions  from  the  hot  bands  of  CO2 
(1/2)  while  detector  31  measured  the  cold  band  radiance.  In  Figs.  1-3  the 
measurements  and  SHARC  predictions  are  compared.  The  measured  radiances 
in  all  three  bandpasses  disagree  with  the  predictions  by  up  to  a  factor  of  1.8. 

Presuming  first  that  the  sensor  calibration  is  correct,  we  estimated  the 
temperature  shift  necessary  to  bring  the  SHARC  predictions  into  agreement 
with  the  data.  This  was  done  using  the  approximation,  valid  for  optically 
thick  radiation,  that  the  ratio  of  the  SHARC  and  observed  radiances  can  be 
represented  by  the  ratio  of  blackbody  intensities  evaluated  at  the  MSIS-90  and 
the  actual  temperatures.  We  calculated  the  temperature  shift  as  a  function  of 
tangent  height  for  each  of  the  three  bandpasses,  from  which  we  calculated  an 
average  shift.  The  temperature  shift  thus  estimated  was  not  uniform,  and  is 
approximated  by  the  following  equations: 
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aT  =  0  if  h  2:  95  km 

=  20  if  h  <  80  km  (1) 

=  20  -  1.33(h  -  80)  if  80  <  h  <  95 


where  h  represents  the  tangent  height. 

The  SHARC  predictions  for  the  CO2  radiances  using  the  new,  shifted 
temperature  profile  are  plotted  in  Figs.  1-3.  The  agreement  with  the  ELC-1 
data  is  much  better  than  with  the  MSIS-90  profile  and  we  conclude  that  the 
modified  temperature  profile  may  a  be  good  approximation  for  the  temperature 
during  the  ELC-1  flight.  The  modified  temperature  profile  is  shown  in  Fig.  4 
as  ELC  (temp).  For  comparison,  we  also  plot  the  unshifted  temperature 
profile  of  MSIS-90. 
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Figure  1.  Comparison  of  CO2  Radiance  Profile  From  ELC-1  Detector 
9  with  SHARC  Predictions.  SHARC  (MSIS)  Assumes  the 
MSIS-90  Temperature  Profile  While  SH>A^C  (temp)  Assumes 
the  Modified  Temperature  Profile  ELC  (temp)  in  Figure  4. 
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Figure  2.  Radiance  Profiles  as  in  Figure  1,  Except  for  Detector  31. 
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Radiance  Profiles  as  in  Figure  1,  Except  for  Detector  7. 
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MSIS  ANO  ELC  TE»«>EflATURE  PROFILES 


Figure  4.  Comparison  of  Derived  Temperature  Profiles  With  the 
MSIS-90  Profile.  ELC  (temp)  Was  Chosen  to  Match  the 
SHARC  Calculation  to  the  Measured  CO2  Radiances.  ELC 
(cal)  Was  Chosen  to  Match  the  Calculation  to  the  Measured 
Radiances  Multiplied  by  1.8. 

An  alternate  explanation  for  the  discrepancy  between  the  predicted  and 
observed  CO2  radiances  would  be  a  calibration  error  in  the  ELC  data 
combined  with  an  error  in  the  kinetic  or  CO2  vibrational  temperature  for 
altitudes  above  80  km.  In  Fig.  4  we  present  the  kinetic  temperature  profile 
needed  to  match  the  SHARC  radiances  predictions  to  "corrected"  data,  if  we 
assume  that  the  original  ELC  data  are  too  low  by  a  factor  of  1.8  for  all 
detectors.  This  profile  is  labelled  as  ELC  (cal).  In  Figs.  5-7,  we  compare 
the  predictions  for  the  ELC  (cal)  temperature  with  the  recalibrated  data.  The 
comparisons  are  as  good  as  those  in  Figs.  1-3. 
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Figure  5.  Comparison  of  1.8x  the  CO2  Radiance  Profile  From  ELC-l 
Detector  9  With  the  SHARC  Prediction  Based  on  the 
Temperature  Profile  ELC  (cal)  in  Figure  4. 
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Figure  6.  Radiance  Profiles  as  in  Figure  5,  Except  for  Detector  31. 
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Radiance  Profiles  as  in  Figure  5,  Except  for  Detector  7. 


The  accuracy  of  the  calibration  of  the  similar  ZIP  instrument  has  been 
discussed  by  Murdock  and  Price.P)  Based  on  their  discussion,  we  believe  that 
the  proposed  calibration  correction  factor  of  1.8  should  be  an  upper  limit. 
We  also  believe  that  the  alternative  possibility  of  a  temperature  shift  of  up  to 
200  C  compared  to  MSIS-90  would  represent  an  upper  limit  to  the  shift, 
given  the  typical  accuracy  of  MSIS-90.  Therefore,  the  two  temperature  profiles 
in  Fig.  4  are  most  likely  the  upper  and  lower  bounds  for  the  true  ELC 
temperature  profile.  These  two  temperature  profiles  are  used  in  the 
subsequent  data  analysis  to  provide  two  separate  estimates  for  the  water  vapor 
profile. 
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3.  H20  data  reduction 


In  this  section  we  discuss  the  procedure  by  which  the  raw  water  vapor 
data  were  reduced.  Pure  rotational  emission  from  H2O  was  measured  by 
detectors  27  and  28  in  the  23  iim  -  29  region.  We  selected  two  sets  of 
data,  one  each  from  the  sets  of  horizontal  and  vertical  scans  mentioned 
earlier.  The  original  data  were  tabulated  in  time  intervals  of  0.003  sec.  In 
order  to  reduce  the  data  points  to  a  reasonable  number,  as  well  as  to  smooth 
out  some  noise,  the  data  were  averaged  over  three  time  steps.  The  maximum 
change  in  tangent  height  in  this  time  interval  is  around  1  km. 

The  next  step  in  the  analysis  was  to  subtract  off  the  background 
radiance  arising  from  off-axis  telescope  leakage.  This  subtraction  affects  only 
the  highest-altitude  water  data  (above  around  80  km).  Assuming  that  above 

110  km  the  signals  in  the  H2O  detectors  are  entirely  due  to  the  leakage 
background,  which  decreases  exponentially  with  tangent  height,  the  magnitudes 

and  scale  heights  for  the  background  were  established  for  each  detector.  The 
data  for  detectors  27  and  28  after  the  background  subtraction  are  shown  in 
Figs.  8-9. 

With  the  background  subtracted,  the  four  sets  of  data  (two  each  from 
the  two  detectors)  were  averaged.  The  results,  which  are  shown  in  Fig.  10 

for  the  limiting  calibration  scale  factors  of  1.0  and  1.8,  were  used  in  all 
subsequent  analysis.  The  profile  is  very  smooth,  since  there  were  sufficient 
data  at  all  tangent  heights  of  interest  and  the  signal-to-noise  ratio  is  high. 

If  the  H2O  radiation  were  optically  thin,  we  could  analyze  the  Fig.  10 

profiles  using  linear  techniques  such  as  Abel  inversion(6)  or  a  simple 
onion-peel  inversion.  However,  as  discussed  in  the  next  section,  H2O  is 
optically  thick  below  85  km,  complicating  the  inversion  process. 
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ELC-1  Water  Vapor  Radiance  Measured  by  Detector  27, 
After  Background  Subtraction. 


CHANNEL  28  H2O  RADIANCE 


CD 


LU 

X 


2 

LU 

CD 


RADIANCE  (W/cm^-sr) 

ELC-1  Water  Vapor  Radiance  Measured  by  Detector  28, 
After  Background  Subtraction. 


-  10  - 


ELC  H2O  (23  4m  -  29  4m) 


tr) 


CD 

z 

<£ 


RADIANCE  (W/cm^-sr) 


Figure  10.  Average  ELC-1  Water  Vapor  Radiance  With  and  Without 
the  Scale  Factor  of  1.8. 
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4.  WATER  VAPOR  PROFILE 


In  order  to  determine  the  opacity  of  H2O  emissions,  we  ran  SHARC 
for  the  ELC'l  temperatures  and  the  standard  SHARC  water  vapor  density 
profile.  At  altitudes  above  85  km,  the  ratio  of  the  full-limb  radiance  to  the 
half-limb  radiance  was  approximately  2,  indicating  optically  thm  emissions. 

Below  85  km,  the  ratio  of  full-limb  to  half-limb  radiances  decreased 

monotonically  to  almost  unity  at  50  km,  indicating  strong  self-absorption. 

4.1  Inversion  Procedures  for  Optically  Thick  Emission 

The  problem  of  inverting  an  optically  thick  limb  emission  profile  is 

essentially  a  non-linear  fitting  problem  and  may  be  solved  by  the  usual 

Newton-Raphson  method.(7)  The  atmosphere  is  approximated  with  discrete 
layers,  and  a  Taylor  expansion  of  the  limb  radiance  is  made  about  an 

initial-guess  profile  and  truncated  after  the  linear  terms.  The  resulting  linear 
equations  are  solved  and  the  process  is  iterated  until  convergence  is  achieved. 
To  reduce  the  dynamic  range  of  the  variables,  and  ensure  that  errors  are 
expressed  in  proportional  rather  than  in  absolute  terms,  we  may  express  the 
number  density  and  limb  radiance  logarithmically.  Then  the  linear  equations 
are  given  by 

In(lj)  =  In(Ioj)  +  ^(dln(lj)/dln(Nj))iln(Ni)  (2) 

i 

Here  j  is  the  index  for  the  line-of-sight  (LOS),  Ij  is  the  limb  radiance  value 
to  be  fitted,  loj  is  the  limb  radiance  calculated  using  trial  densities  Noi  in 

each  layer  i,  Nj  is  the  density  in  layer  i,  and  Aln(Ni)  =  ln(Nj)-ln(Noi).  The 
species  density  may  refer  to  either  ground  state  or  excited  state  number 
density,  if  it  is  assumed  that  the  relationship  between  them  is  governed  by  a 
known  kinetic,  vibrational,  or  other  temperature  in  each  layer. 

The  partial  derivatives  in  Eq.  (2)  correspond  to  the  usual  weighting 
functions  and  may  be  computed  by  incrementing  the  layer  densities  that  are 
input  to  the  radiation  code.  The  Nj  are  obtained  by  solution  of  the 


simultaneous  equations,  which  can  be  performed  by  successive  substitution 
starting  from  the  top  layer.  To  achieve  convergence,  successive  iterations  can 
be  performed  using  the  Nj  as  the  new  trial  values  NqI- 

In  the  current  SHARC  code,  it  is  not  possible  to  increment  the  layer 
densities  without  generating  a  new  population  file,  which  requires  re-executing 
the  chemistry  and  radiation  transport  (CHEMKIN  and  NEMESIS)  routines. 
Determination  of  the  weighting  functions  in  this  manner  would  be  very 
cumbersome  and  time-consuming.  However,  an  external  or  internal  routine 
could  be  developed  to  automatically  update  the  population  file  and  execute  the 
LOS  radiance  computation  alone. 

For  the  present  application,  we  have  chosen  a  simpler  inversion 
approach  that  minimizes  the  computational  effort  and  uses  only  the  existing 
SHARC  code.  Accurate  forward  calculations  are  performed  for  each  trial 
profile,  while  the  partial  derivatives  that  drive  successive  iterations  are 
approximated  using  the  crude  assumption  that  the  distribution  of  the  species 
column  density  along  the  LOS  is  immaterial.  This  would  be  the  case  if  all 
temperatures  (kinetic,  vibrational,  etc.)  were  uniform  along  the  LOS.  We  can 
then  replace  Eq.  (2)  by 

In(Ij)  =  ln(I„j)  +  (dln(Ij)/dln(Nj))Aln(Nj)  (3) 

where  Nj  is  the  column  density  along  the  j’th  LOS  and  Aln(Nj)  s 
ln(Nj)-ln(Noj),  Ngj  being  the  trial  column  density.  From  Eq.  (3),  the  Nj  are 
given  by 

Nj  =  Noj(Ij/loj)>/“j  (4) 

where  aj  s  dln(Ij)/dln(Nj).  The  Nj  can  be  inverted  to  determine  the  layer 
densities  Nj  using  a  standard  linear  inversion  method.  For  optically  thin 
radiation,  a  =  1  and  this  method  reduces  to  a  linear  inversion  of  the  column 
radiance.  For  optically  thick  radiation,  a  <  1,  which  implies  that  for  a  given 
fractional  change  in  species  density  the  associated  fractional  change  in  limb 
radiance  is  smaller. 

The  aj  can  be  determined  from  only  two  sets  of  SHARC  calculations 
with  different  column  density  profiles,  such  as  full-limb  calculations  using 
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Table  2.  Values  of  a;  =  dln(Ij)/dln(Nj),  Where  Nj 

is  Water  Vapor  Colunm  Density  Along  the  LOS. 


Tangent  Ht 
(km) 

Half-Limb  Method 
(Standard  Profile) 

Full-Limb  Method 
(Standard  and  First 
Iteration  Profiles) 

>  90 

1.00 

1.00 

85 

0.97 

0.97 

80 

0.80 

0.89 

75 

0.53 

0.68 

70 

0.36 

0.44 

65 

0.25 

0.34 

60 

0.19 

0.37 

55 

0.18 

0.38 

50 

0.23 

0.46 

successive  density  profile  iterations  or  full-limb  and  half-limb  calculations  using 
the  same  density  profile.  We  used  values  derived  via  the  half-limb  method, 

which  requires  only  a  single  CHEMKIN  and  NEMESIS  calculation,  to 
determine  the  first-iteration  profile,  and  used  this  profile  and  the  original  guess 
to  derive  refined  a  values  that  were  used  in  subsequent  iterations.  The  a 

values  from  both  methods  are  given  in  Table  2. 

Due  to  the  uniform-LOS-temperature  assumption  in  the  above  inversion 
scheme,  the  procedure  may  not  be  fully  convergent  for  realistic  atmospheres. 
However,  the  error  can  be  minimized  by  a  judicious  definition  of  N.  For 

example,  if  the  emission  is  due  to  a  vibrationally  excited  state,  choosing  N  as 
that  state  rather  than  the  ground  state  removes  the  first-order  error  in  the 
optically  thin  limit.  For  the  present  case  of  rotationally  excited  states,  we 

define  N  as  an  "effective"  excited  state  species  column  density, 

Nj  =  2Ni‘ij«’‘P(-Ea/kTi)  (5) 

i 

Ea  is  the  average  excited  state  energy  (approximately  1100  cm*l  for  the  strong 
water  rotation  lines  in  the  ELC-1  bandpass),  Tj  is  the  temperature  in  layer  i, 
Ijj  is  the  length  of  the  LOS  through  layer  i,  and  Nj  is  the  water  vapor 
number  density  in  layer  i. 
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In  conclusion,  the  inversion  procedure  given  by  Eqs.  (3)-(5)  should 
provide  accurate  solutions  in  optically  thin  conditions  (a  =  1)  as  well  as  with 
a  uniform  temperature  profile  in  optically  thick  conditions.  As  shown  in  the 
following  subsection,  good  results  have  been  obtained  in  the  current 
investigation  of  water  vapor  rotational  emission.  However,  complete 
convergence  under  general  conditions  may  require  a  fuller  treatment  such  as 
given  by  Eq.  (2).  The  development  of  an  inversion  algorithm  along  these 
lines  utilizing  the  SHARC  code  would  be  very  desirable  for  future  work. 

4.2  Application  to  ELC-1  Water  Vapor  Data 

The  inversion  procedure  based  on  Eqs.  (3)-(5)  was  implemented  using 
the  ELC-1  limb  radiance  profiles  in  Fig.  10.  With  one  or  two  iterations, 
convergence  of  the  observed  and  calculated  radiances  was  obtained  to  within 
10%  or  better  at  all  tangent  heights  from  55  km  to  85  km.  The  results  are 
shown  in  Figs.  11-12.  In  Fig.  11,  the  comparison  is  to  the  predicted  radiance 
for  the  temperature  profile  labelled  ELC  (temp)  in  Fig.  4.  In  Fig.  12,  the 
comparison  is  for  the  temperature  profile  labelled  ELC  (cal).  In  both 
comparisons,  the  agreement  is  excellent.  It  is  not  clear  whether  the  residual 
error  is  due  to  a  shortcoming  of  the  inversion  method  itself  or  its  method  of 
implementation.  The  comparisons  between  the  calculated  and  observed 
radiances  in  Eq.  (4)  were  made  at  5  km  tangent  height  increments,  while  the 
SHARC  code  uses  2-km  layering.  Given  the  very  small  scale  height  of  water 
vapor  emissions,  finer  altitude  grids  may  ultimately  be  required. 

In  Fig.  13,  we  compare  the  two  ELC-1  water  vapor  density  profiles  (for 
the  two  temperature  profiles)  to  the  SHARC  standard  profile.  At  the  lower 
altitudes,  the  true  density  profile  should  lie  between  the  two  profiles  shown. 
At  higher  altitudes,  the  profiles  coalesce.  However,  there  the  effect  of  tangent 
height  uncertainty  must  be  included  in  an  overall  error  estimate.  From  the 
data  in  Fig.  9,  a  2  km  tangent  height  spread  is  associated  with  a  radiance 
spread  of  nearly  a  factor  of  2,  and,  at  these  altitudes,  a  comparable 
uncertainty  in  water  density.  Considering  the  temperature,  tangent  height  and 
inversion  uncertainties,  reasonable  bounds  for  the  true  water  vapor  profile  in 
ELC  are  obtained  by  averaging  the  two  profiles  and  taking  factor-of-two  error 
bars  at  all  altitudes.  The  result,  expressed  as  a  mixing  ratio,  is  plotted  as  the 
shaded  area  in  Figs.  14  and  15. 
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Figure  11.  Comparison  of  Measured  Water  Vapor  Radiance  Profile 
With  the  SHARC  Simulation  Using  the  Density  Profile 
ELC  (temp)  in  Figure  13. 
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Figure  12.  Comparison  of  Measured  Water  Vapor  Radiance  Profile, 
Including  the  Scale  Factor  of  1.8,  With  the  SHARC 
Simulation  Using  the  Density  Profile  ELC  (cal)  in  Figure 
13. 
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Figure  13.  Comparison  of  the  ELC-1  Water  Vapor  Profiles  and  the 
SHARC  "Standard"  Profile. 
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Figure  14.  Comparison  of  the  ELC-1  Water  Vapor  Mixing  Ratio 
(Shaded  Area)  to  Space-Based  IR  Absorption  Measurements 
Made  During  the  Fall/Winter  Months. 
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Figure  15.  Comparison  of  the  ELC-1  Water  Vapor  Mixing  Ratio 
(Shaded  Area)  to  Ground-Based  Microwave  Measurements 
Made  During  the  Fall/Winter  Months. 
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5.  CONCLUSIONS 


In  Fig.  14  the  ELC-1  water  vapor  profile  is  compared  with 

measurements  from  three  ground-based  microwave  experiments  at  mid-latitudes. 
These  determinations  are  based  on  analysis  of  the  22.2  GHz  emission  line 
shape,  and  are  limited  to  altitudes  below  70-80  km.(^)  The  microwave 
measurements  were  made  at  three  different  sites:  Haystack  Observatory  in 
Westford,  MA,(9)  JPL,(1'^)  and  Penn.  State.(ll)  Measurements  have  been 

made  at  various  times  of  the  year  and  a  significant  seasonal  variation  has 
been  observed.  The  water  vapor  mixing  ratio  has  been  seen  to  vary  by  an 
order  of  magnitude  at  the  higher  altitudes  and  by  about  a  factor  of  two  at 

the  lower  altitudes,  and  is  greatest  during  the  late  summer  months. 

In  Fig.  15  we  compare  the  ELC-1  result  to  two  other  measurements 
from  space  platforms.  These  were  made  during  the  first  Spacelab  mission  on 
board  the  Space  Shuttle,  in  the  winter  of  1983(^2)  and  from  a  rocket 

launched  from  Cape  Perry  on  December  6,  1974.(13)  in  both  experiments, 

the  solar  absorption  spectrum  of  the  2.7  nm  band  of  water  vapor  was 

measured  during  sunrise  or  sunset. 

Figures  14  and  15  show  that  the  water  vapor  profile  derived  from  the 
ELC-1  data  is  consistent  with  previous  mid-latitude  measurements  made  during 
fall  or  early  winter  months.  The  uncertainty  in  an  individual  profile  is  about 
a  factor  of  two  and  the  reader  is  referred  to  the  original  reportings  for  the 

exact  uncertainties.  The  range  of  these  measurements,  while  substantial,  is 
much  smaller  than  for  measurements  that  are  not  selected  according  to  season. 

The  results  from  ELC-1  demonstrate  that  optically  thick  rotational 
emission  can  be  analyzed  to  yield  the  altitude  profile  for  water  vapor  in  the 
mesosphere  with  an  accuracy  comparable  to  microwave  and  IR  absorption 
methods.  Global  data  from  CIRRIS-1  can  be  analyzed  in  a  similar  manner 

and  can  be  expected  to  provide  valuable  information  for  understanding  and 
modeling  water  vapor’s  variability. 
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